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ABSTRACT 

The primary purpose of this contribution is to develop a novel framework for generalized 
robust design of tuned mass damper (TMD) systems as passive vibration controllers for 
uncertain structures. This versatile strategy is intended to be free of any restriction on the 
structure-TMD system configuration, the performance criterion, and the number of 
uncertain parameters. The main idea pursued is to adopt methods and concepts from the 
robust control literature, including: (1) the linear fractional transformation (LFT) formulation 
pertaining to the structured singular value (;<) framework; (2) the concept of weighted multi- 
input multi-output (M1MO) norms for characterizing performance; and (3) a worst-case 
performance assessment method to avoid the unacceptable computation burden involved 
with exhaustive search or Monte Carlo methods in the presence of multiple uncertainties. 
Based on these, the robust design framework is organized into four steps: (1) modeling and 
casting the overall dynamics into the proposed LFT framework that isolates the TMD system 
as the controller, and the uncertainties as a structured perturbation to the nominal dynamics; 
(2) setting up the optimization problem based on generalized indices of nominal 
performance, robustness, and worst-case performance; (3) implementing a genetic algorithm 
(GA) for solution of the optimization problem; and (4) post-processing the results for syste- 
matic visualization, validation, and selection of preferred designs. This strategy has been 
implemented on several illustrative design examples involving a seismically excited multi- 
story building with different combinations of assumptions on the uncertainty, TMD con- 
figuration, excitation scenarios, and performance criteria. The resulting solution sets have 
been studied through various post-processing methods, including visualization of Pareto 
fronts, uncertain frequency response plots, time-domain simulations, and random vibration 
analysis. 

© 2010 Elsevier Ltd. All rights reserved. 



1. Introduction 

Due to advantages of efficacy, reliability, and low cost, TMD systems have been extensively studied for implementation on 
various structures, such as buildings [1,2], bridges [3], and tension leg platforms [4], Some possible TMD realization schemes 
for controlling horizontal motion of buildings (including spring-damper, rubber bearings, and various pendulum schemes) 
have been discussed in [5]. Sometimes innovative physical realizations are also possible, e.g. see the roof-garden TMD design 



* Corresponding author. 
E-mail addresses: amohtat@ut.ac.ir, amohtat@cim.mcgill.ca (A. Mohtat), edehghan@ut.ac.ir, edehghan_niri@yahoo.com (E. Dehghan-Niri). 
1 Current address: Department of Mechanical Engineering and Centre for Intelligent Machines, McGill University, 817 Sherbrooke St. West, Montreal, 
Quebec, Canada H3A 2K6. 

0022-460X/S-see front matter © 2010 Elsevier Ltd. All rights reserved, 
doi: 1 0.1 01 6/j.jsv.201 0.09.007 



A. Mohtat, E. Dehghan-Niri / journal of Sound and Vibration 330 (2011) 902-922 



Nomenclature 




space realization: 


Important symbols 




x = Ax+Bv 
y = Cx+Dv 


d 
e 


disturbance input vector, e.g. ground accelera- 
tion and wind loading 

normalized control-target (error) output vec- 
tor 


lower and upper LFTs defined as (provided that 

the matrix inverses exist): 
r e (M,N) ^(M,N) = M„ +M 12 N(I-M22Nr 1 M 2 i 
^ U (M,Q) ^ U (M,Q) = M 22 + M 2 i Qfl-IVl! ! QF 1 M 12 


F c 


control force input to the plant exerted from 
the TMD system 




where NsC 2 ^ 2 , Q. e C" 1 xPl and Me 

c (p,+p 2 )x Wl +<) 2 ) is cons i stent iy partitioned 


JNP. JR 


and Jwcp indices of nominal performance, 








robustness, and worst-case performance 


List of 


acronyms 


K(s) 


TMD system as the controller 






P(s) 


plant dynamics, either 2-in/2-out ports or 3- 


DMF 


dynamic magnification factor 


V 


in/3-out ports 

normalized disturbance input vector 


FRF 
GA 


frequency response function 
genetic algorithm 


W;(s), W (s) input/output weighting filters 


LFT 


linear fractional transformation 


y 

z 


feedback output vector from the plant (input 

to the controller) 

control-target (error) output vector 


LMI 

MIMO 

SDOF 


linear matrix inequality 
multi-input multi-output 
single degree of freedom 


,5 

A 


uncertain real scalar element, - 1 < S < 1 
structured uncertainty perturbation block 


SISO 
SSV 


single-input single-output 
structured singular value (/j.) 






TITO 


two-port input, two-port output 


Some mathematical notations and definitions 


TMD 


(also: STMD, MTMD, UMTMD, IMTMD) tuned 








mass damper, also: single TMD, multiple 




A 


B 


shorthand for system transfer matrix 
C(sl-A) ] B+D, corresponding to the state- 




TMD (MTMD), uniform MTMD, irregular 
MTMD 


C 


D 



in the aforementioned reference. From a theoretical point of view, on the other hand, the STMD can be regarded as a narrow- 
band filter tuned to the first structural frequency. This leads to vulnerability to mistiming effects [6] and degradation of the 
STMD efficacy in mitigating structural response against excitations with a wide spectrum of frequency contents, like seismic 
loads [2]. Moreover, a single massive TMD is rather difficult to install and maintain [7]; and, when subjected to dynamic 
impulsive loads (like earthquake excitations) filtered through the structure, does not reach its full energy-absorbing potential 
as quickly as required [2]. These factors have motivated the application of systems with several TMD units covering a 
frequency range and/or physically distributed in space. In this paper, the term "TMD system" refers to any such configuration 
of a number of TMD units installed on a generic structure. Examples of these configurations include the dual TMD [8] also 
referred to as bi-TMD [9]; double TMD [10]; double-layered TMD [10]; interconnected multiple TMD [11] (interconnection of 
TMDs in series); and, non-interconnected MTMD (NI-MTMD) [11] or simply referred to as MTMD which is the most broadly 
used configuration. Different assumptions on the parameters have been considered for design of TMD systems. For instance, 
in the MTMD case, equal masses and damping ratios along with uniform distribution of frequency ratios can be assumed for 
the TMD units to simplify the design; resulting in a uniform MTMD (UMTMD) [12] compared to the irregular MTMD 
(IMTMD) [13] where frequency and damping ratios of the TMD units constitute individual design variables. As previously 
pointed out, some other investigations study physical distribution of TMD devices on structures, e.g. references [1,2] study 
optimal placement of TMDs on building floors. Despite this extensive research, there is almost no study that offers a 
convenient method and formulation to accommodate any arbitrary combination of structure-TMD system within a unified 
approach. Other limitations of the investigations available in the literature are discussed in what follows. 

In practice, structures and their dynamic characteristics are subjected to uncertainties in the mass, damping and stiffness 
originating from a variety of sources, such as inadequate modeling of the boundary conditions at the structural joints, effects of 
nonstructural elements, degradation due to aging, and fluctuations in structural mass, as well as uncertainties in the member 
capacities, yield strengths, etc. Furthermore, the dynamic characteristics of structures change under earthquake or wind 
excitations. Some structures have uncertain nonlinear properties even in the small amplitude range due to the contribution of 
secondary members [15,16]. Therefore, in practical design problems, robustness of the TMD system to uncertainties must be 
treated as a design concern to avoid undesired performance degradation and mistuning effects. Nevertheless, in most of the 
previous works, optimization of the TMD systems was based on effectiveness; and, robustness was either not considered [17- 
21 ] or qualitatively checked and commented on, only after the design [7,12,22-25]. See the brief discussion in [16] with further 
details on these studies and the qualitative concepts regarding robustness. However, quite recently more attention has been 
dedicated to the issue of robustness in TMD systems. Hoang and Warnitchai [26] proposed a robustness index defined as an 
average performance of several samples of the uncertain system with the uncertain parameters following a specified 
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probability distribution. Their procedure, although rather extendable to situations more general than MTMDs for SDOF 
structures (the numerical example considered in that paper), is limited to mean square structural response reduction against 
stationary stochastic excitation due to the Lyapunov-based formulation of the performance function and its gradient. In another 
case in point, Li and Ni [27] proposed a rigorous gradient-based method that has the ability to design a robustly effective 
IMTMD considering a desired level of estimation error in structural parameters. However, this procedure was intended only for 
minimization of the maximum DMF of structural response (SISO case of H„_. optimization) and was presented only for a SDOF 
structure with MTMD. In another relevant research, Dehghan-Niri et al. [16,28] proposed the novel idea of utilizing an index of 
robustness along with effectiveness within a two-objective evolutionary optimization and demonstrated the efficacy of this 
methodology for design of uniform and irregular MTMD systems for SDOF main structures with uncertain natural frequency. 
Similar two-objective GA design methodologies have been proposed in [9,29]. In Ref. [29], Marano et al. proposed the 
minimization of the mean and standard deviation of the structural displacement covariance computed through sampling as the 
two objectives. This method aimed only at vibration mitigation of SDOF main systems with STMD against random base 
acceleration modeled by a stationary filtered white noise process. In Ref. [9], Ok et al. demonstrated the application of GA with 
nominal and worst-case robust performance indices as the two objectives for design of bi-tuned mass dampers attached to the 
SDOF main system with uncertain stiffness. The main purpose of the aforementioned investigation was to develop practical 
optimal design formulae for bi-TMDs within methodology and formulations that maintain simplicity. 

In conclusion, each of the available methodologies for design of TMD systems targets a specific purpose and has its own 
limitations and restrictive assumptions. These limitations can be categorized as: 

1 ) Restrictions on the structure and TMD system configuration, e.g. analytical formulations assuming a SDOF structure 
with a MTMD. 

2) Restriction on design or performance criteria, e.g. a specific SISO norm from a specific excitation to a specific target output. 

3) Considering only a simplistic single-parameter uncertainty, or treating the multiple-parameter uncertainty by 
inefficient Monte Carlo-like (sampling) methods. There exist only few studies that use more computationally efficient 
asymptotic expansion methods for reliability integrals (e.g. see Ref. [15]). Furthermore, these investigations target only 
aleatory (stochastic) uncertainties (see Section 4 for a categorization of uncertainty types), and are constrained by 
conditions to maintain the required analytical formulations tractable. 

The present paper, on the other hand, is intended to overcome all these limitations. The first one is handled by the lower 
LFT formulation introduced in Section 2; the second is resolved via the concept of weighted MIMO norms from any group 
of disturbance v to error outputs e; and, the third is tackled through explicit representation of the multiple uncertainties as 
a structured perturbation block within the rigorous LFT formulation of Section 4.2. This section exploits the SSV (/i) 
framework [14] and the LFT-based worst-case performance assessment methods [30] in the robust control literature. To 
facilitate the design of the TMD systems in such a generalized sense, all the concepts, formulations and procedures are 
collected into the four-step framework shown in Fig. 1. 

It should be pointed out that the organization of the paper is not based on the four steps of the framework, but on the 
logical flow of the concepts. The complete development of the LFT formulation required in Step I is postponed to Section 4.2 
after its initial development in Section 2. The procedures associated with Step II are discussed in Sections 3 and 4, while Step 



Step I 






Step II 






Modeling & Casting the Dynamics 
into the LFT Formulation 

• Modeling the structure as the plant 

• Selection of uncertain parameters 
and pulling out their 8's 

• Modeling the TMD system as the 
controller 




Optimization Problem Setup 

• Choosing design variables and 
constraints 

• Selection of in/out weighting filters 

• Choosing the objectives and 
writing subroutines for computing 
them 












V 




Post-processing of Results 

• Visualization (uncertain FRF plots, 
Pareto fronts, etc.) 

• Validation (time-domain 
simulations, random vibration 
analysis, etc.) 

• Selection of preferred designs 












GA Implementation 

• Algorithm type 

• GA tuning parameters 

• Convergence monitoring 










Step III 





Step IV 
Fig. 1. Flowchart of the proposed framework for robust design of TMD systems. 
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IV is practically illustrated throughout the examples worked out in Section 5. GA (Step III), on the other hand, is considered a 
standard optimization tool; hence, there is no specific section dedicated to it. Of course, some important relevant concepts 
are briefly discussed in Sections 3-5, because the proposed framework greatly owes its simplicity and flexibility to GA. 



2. LFT formulation for the structure-TMD system 

In this section, we shall show that the dynamics of the structure with the TMD system can be cast as the LFT 
interconnection of a TITO plant and a feedback controller, as shown in Fig. 2. 

Considering disturbance d and control F c as the two groups of inputs, and error output z and feedback y as the groups of 
outputs, the TITO model for the structure referred to in Fig. 2 is 



P(,) = 



A 


B, 


B F 


C 


D. 


D, 


c,. 


D, rf 


n yF 



while the TMD system model is 



K(S): 



\A K 


B,l 


c K 


D *_ 



(1) 



(2) 



and the lower LFT [14] interconnection describing the mapping from the input group of signals d to the output group z would be 

T zd (s) = ^-,(P.K) = Pzd + P zf K(I-P yf K) J P yd (3-1) 

where 



p„(*) = 



A 


Bj 


C. 

"A 


D '_ 
B/ 


c, 


D '. 



P=rW = 



, VJ S ) : 



A 


B,l 


A 


B/ 


c r 


V 



are the partitions of the transfer matrix P, consistent with the input/output grouping; or, equivalently: 



T ;(/ (,)=^(P,K): 



A-CL 


B ci ] 


^CL 


D CI 



(3-II) 



where A C l, B C l, C C l, and D C l are the state-space matrices of the closed-loop system, defined as 

A+B F R 'DjfCy B F R _1 C K 



Bi^R C v 



A«: + Bj<-R C yf Cj<- 



B d + B f R D K C 



K^yd 



B/fR C yt j 



: CL = Q + C^DKRCy QpR-'Q] 



D CL = D zd + D zF D K R D 



•yd 



with R -. 



CypDn, R ■ 



-DKCyF. 







Structure 

as the 

TITO Plant P 






















F c 






y 




TMD System 

as the FB 
Controller K 















Fig. 2. LFT interconnection of the structural plant and the TMD controller. 
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This formulation isolates the dynamics of the TMD system to be designed from that of the structural plant, facilitating the 
performance evaluation of designs. It offers a systematic convenient way for definition of weighted M1MO norms as a means of 
characterizing the performance of the controlled structural response from disturbance inputs to target quantities (see Section 
3). The LFT machinery will be further exploited in Section 4.2 to pull out the parametric uncertainties embedded in the system. 



2. 1 . TITO formulation for the structure 



A generic scenario of an N-story shear building subjected to base acceleration and force excitation, as shown in Fig. 3, is 
utilized to illustrate the required formulation for the structure. 
The matrix equation of motion for the structure is as follows: 



(4) 







M s u + C 5 u + K s u = 


-M s r g a g + F d + F ctr i 




Zi = W u u, 


z 2 =ii + a g 


where 






M s =diag(m 1 , m 2 , . 


. , m N ) is the structural mass matrix: 








~k 1 +k 2 -k 2 









-k 2 k 2 +k 3 


-*3 




Ks = 




— k N _] k N _i+kn —kn 









-k N k N 



is the structural stiffness matrix. 

2^Mia>2 ivr j — hi — k s is the proportional damping matrix (any other dissipation matrix can be also considered 



s OJi + C02 S (Jj\+0] 2 ' 



without additional effort, if required) defined based on the damping ratio f s and the first two natural frequencies toi, co 2 of 
the structure. Note that the dissipation matrix assuming that damping forces depend only on generalized velocities is not 
the only linear model of vibration damping (see Ref. [31] for a detailed discussion). However, for simplification of the 
illustration, a dissipation matrix damping is assumed for the shear building model (this does not affect the generality of the 
framework, because the performance measures are formulated using the plant P without any assumptions rather than 
linearity). 

r g =[l ...1] T is the influence vector of the ground acceleration a g . 

Fd . Fctri are the disturbance forces (e.g. from wind load) and control forces (exerted by the TMD system) to the floors of 
the building. 

1 01 



W„ 



-1 1 



is the weighting matrix utilized to extract the inter-story drifts (floors displacement relative to 



-1 1. 

each other) z 1 = [u 1 , u 2 — u-[ u N -u N _ 1 ] 1 out of the floor displacements u relative to the base 

and 

z 2 is the vector of floors absolute accelerations. 




Fig. 3. JV-story shear building subjected to base acceleration and force disturbances, equipped with "MTMD on i'th floor" (a) and "STMD on every floor" 
(b) systems. 
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Now, considering the following definition for the state vector and the input-output groups of the TITO model: 
states: x= ■! . )•; inputs: d={ ° }, F c = F ctr i; outputs: z= < }, y = z 2 



the following definition for the state-space matrices previously used in Eq. (1) would result: 

OnxN InxN 



-Mj'Ks 


-*£•'£ 


w u 


OnxN 


M^K S 


-M^Q 


-M^K* 


-Mr'C 



B,, 



D 



'yd = 



N xl 

-'g 

"nxI 
"nxI 



OnxN 
M s _1 


, B F = 


" 


JnxN 




®NxN 

M7 1 


, D zF = 


OnxN 

M7 1 



m; 



D yF = [MJ 



(5) 



which completely determine the system P from disturbance d (base acceleration and force excitation) and control F c input 
(TMD system forces exerted to the floors) to error output z (inter-story drifts and floors absolute acceleration) and 
feedback output y (absolute acceleration of building floors). 



2.2. TMD system as the feedback controller 

The equation of motion for a single TMD unit in terms of u T , its displacement relative to the ground, can be easily 
derived as 



m T (UT + a g ) = mj(UT-Ui + ai) = -k T (UT-Ui)-c T (UT-Ui) 



(6) 



where m T , k T , and c 7 are the TMD mass, stiffness, and damping constant; u, is the displacement of the TMD support (the fth 
floor to which the TMD is attached) relative to the ground; and, a, is the absolute acceleration of the ith floor (a,- = a g + u{). 
The right hand side of this equation is the sum of the spring and dashpot force to the TMD mass, the opposite of which is 
exerted to the TMD support (ith floor) and is denoted by F,. Therefore, the TMD unit can be regarded as a system block from 
the support acceleration input (a,) to the force output (F,), described by 



x T =A r x r + B r a f ; x T =[u T —Uj u T —Uj 

F; = CjXj 



(7) 



where 





-k T /m T 



1 

-c T /m T 



Br 



-1 



[k T c T ] 



Accordingly, each of the TMD systems configurations shown in Fig. 3 (as well as any other possible) can be viewed as an 
interconnection of single TMD blocks as shown in Fig. 4. From the block-diagram representation of Fig. 4, it can be easily 
concluded that the state-space matrices of the TMD system as the feedback controller K in Eq. (2), for the n-MTMD 
installed on the ith floor case can be written as 



[A,d; 



.blkdiagi A T1 ,A T2 A Tn ) 



(a) 












(b) 
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Fig. 4. (a) MTMD interconnection and (b) "STMD on every floor" interconnection. 
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[B K ]; 



B72 

Bt„ 



G a ; Ga = \ga(])h, 



ga(j) ■- 



Vo 



] = * 

otherwise 



Krl 



= G F [C n C T2 ... Cm]; G F = G a r 

[DkInxN = OnxN (8) 

where blkdiag refers to the block diagonal arrangement of matrices. On the other hand, for the "STMD on every floor" 
configuration in Fig. 4b we would have 

[Ajf]2nx2n = blkdiag(A n ,A r2 , 
[B K ] 2nxN = blkdiag(B n ,Br2- 
[ c K] Nx2n = blkdiag(C n ,C T2 ,. 



..Am) 

..B Tn ) 
.C rn ) 



Pit], 



Oiv 



(9) 



Note that n refers to the number of TMD units, while N is the number of building floors; and, these two coincide in the case 
of "STMD on every floor". Any other distribution of TMD units on a building can be also formulated in a similar manner. 
This approach also facilitates the redundancy study of the TMDs (effect of failure of individual TMD units on the overall 
performance). 



3. The nominal performance optimization problem 



In the previous section it was shown that in view of the LFT formulation, we can consider the TMD system as a 
controller for the structure that is intended to improve the closed-loop system response to disturbances. Consequently, a 
natural performance objective for designing this controller would be reducing the closed-loop gain from disturbance 
excitations d to target outputs z. System norms provide a precise measure for this purpose. For a MIMO system, as is the 
case in point, a weighted matrix system norm would be preferred to account for: (1) relative magnitude of different 
disturbance signals composing d; (2) frequency/spectral content of the signals; and (3) relative importance of error outputs 
z in the design. The input-output weighted configuration for this purpose is shown in Fig. 5. 

Accordingly, minimization of the following generalized weighted nominal performance index would be an appropriate 
criterion for tuning the TMD system design variables x to achieve optimal nominal performance: 



miny NP = ||T e „(s)|| = |W (s)T zd (s)W ; (s)|| = |W (s)^(P(s),K(s))W,(s)| 



s.t. design constraints 



(10) 



where ||.|| p refers to the H 2 or H rr . 



system norm [14] (p=2 or 00), while W, and W are the input and output weights 
selected as proper stable transfer matrices of dimensions dim(d) x dim(v) and dim(e) x dim(z), respectively (for the 
problem formulated in Section 2.1: dim(d)=N+l, dim(z)=2N). Formal definitions, existence conditions and computational 
procedures for these H 2 and H, weighted MIMO norms have been discussed in the Appendix. 

As previously pointed out, the input weighting transfer matrix W; enables us to take into account the relative 
magnitude and frequency content of the disturbance signals d consisting of base excitation (earthquake acceleration) and 
force loading (wind-induced forces to the floors), with a general form (assuming dim(v)=N+l) of: 



[W,(s)], 



(N+l)x(N+l) " 



WeqOOlxl 

OnxI 



OlxN 
W ext (S) N 



(11) 



where W eq is the earthquake weighting filter and W ext is the weighting for the external forces to the floors. If the main 
purpose of the TMD system is vibration attenuation against earthquake, as is the case in the examples considered in this 
paper, W ext can be set to zero; then, the all-zero columns can be eliminated to yield the following simplified input weight: 



[W,; 



(N+l)xl • 



:[W eq 0,, 



(12) 
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Fig. 5. Input-output weighted LFT configuration. 
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with dim(v)=l (only base acceleration) and the earthquake weighting simply chosen as either 1 (no frequency content 
information) or 



2Cga>gS+ct>! 

' S 2 +2C„COgS + Ct>2 



We,W = W0 ' ' « d3) 



where co g and f g are characteristic ground frequency and damping ratio, and w is a scaling factor. This is a Kanai-Tajimi 
[32-34] shaping filter that leads the designed TMD system better be trained (during the optimization process) against an 
earthquake loading with a similar spectral content. The effect of this spectral content training will be illustrated in Section 
5.1.2 with details and results. More complicated filters, such as the three stage filter used in [35], that enable a more 
detailed modeling of the ground dominant frequency and soil characteristics can also be similarly incorporated. Moreover, 
details on across-wind-excitation spectral loading can be found in the literature, e.g. [36,37], and can be used to form 
appropriate input weights for design against wind loadings. 

The role of the output weight, on the other hand, is to specify the amount of the designer's emphasis on reducing each of 
the different error outputs z consisting of the inter-story drifts and absolute accelerations of the floors. The general form of 
this weighting, assuming dim(e) =2N, is 



[W (s)] 2 



W drf (s) NxN 



(14) 



with W drf and W acc as the drift and acceleration output filters. But, if the main focus is reduction of inter-story drifts 
(protection against structural damage), as is the case in all the design examples of Section 5, the following simplified 
output matrix weight (after elimination of all-zero rows) is a suitable choice: 

[ W Olkx2N = [ l NxN 0(VxN] (15) 

with dim(e)=N (only inter-story drifts as target controlled quantities); otherwise, if the primary emphasis is reduction of 
floors absolute acceleration (protection against losses to the building contents), a simple choice would be as follows: 

[W 2]nx2N= [°NxN Inxn] (16) 

again with dim(e)=N but with the accelerations selected as target outputs. 

The optimization problem, i.e. Eq. (10), might be solved by different methods. To overcome difficulties and 
disadvantages involved in methods such as exhaustive search (unreasonably low efficiency and high computational cost) 
and gradient-based methods (need for analytical formulation of the gradients, involved algorithmic complexities and 
convergence concerns that seriously restrict the span of treatable problems), genetic algorithm (GA) has been selected. 
Inspired from the concept of natural selection, GA shows a great promise of convergence to a nearly global optimum. GA 
can also easily handle design constraints in Eq. (10): bounds on design variables x are naturally incorporated, while other 
potential constrains on the TMD systems can be treated by penalty functions. Besides searching a descritized hypercube of 
real design variables, GA has also the ability to search a combinatorial space. This is due to the flexibility in the genetic 
coding and the independence to explicit parameterization of the objectives by design variables. This means that except for 
the TMDs specifications, the number, location and configuration of the TMD units can also be simultaneously optimized. Of 
course, exploiting all these possibilities is not the main purpose of this study. Due to these and other many advantages, 
genetic algorithms are well-established in the literature (see Ref. [38] for details on GA operators and implementation 
concepts), and widely used in various engineering problems including optimization of TMD systems [9,16,18,29]. 

4. The robust performance optimization problem 

In the single objective optimization of the previous section, a perfect match between the model and actual structure 
was assumed that if violated might result in a significant degradation of the actual effectiveness of the designed TMD 
system. However, as explained in Section 1, this perfect match cannot be expected in practical situations due to the 
existing uncertainties. Uncertainty is generally categorized as aleatory (due to inherent variation) and epistemic (due to 
lack of exact knowledge). The former is usually modeled by random variables following specific probability distributions, 
while the latter is modeled by parameters each lying in a specific interval [39]. This section is intended to tackle robust 
optimization in the presence of generic multiple epistemic (interval) uncertainties. Section 4.1 deals with the simplest 
form of the uncertain problem assuming only a single real uncertain parameter. This provides a link to the previous 
robustness concepts in the literature and serves as a preamble to the much more complicated case of multiple-parameter 
uncertainty studied in Section 4.2. 

4.1. Single-parameter uncertainty 

The mismatch between natural frequencies of the actual structure and its model has a most significant detuning effect 
on the TMD system effectiveness. Hence, in this subsection, the single-parameter uncertainty is chosen as the following 
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uncertain perturbation to the structural stiffness matrix: 

K sS = K s , nom (1+w K S), -1<<5<1 (17) 

where K s nom is the nominal stiffness matrix, K sS is the perturbed structural stiffness matrix, w K is the uncertainty level, and 
S is the real parametric uncertainty. We shall also modify our notation P for the structural TITO system to P nom for the 
system resulting from K s , n om and P s for that resulting from K sS ; and consequently, a generalized weighted performance 
index for the uncertain structure can be defined as 

J P (<5)=||W F f (P (S ,K)W / || p (18) 

based on which the previously defined nominal performance index / NP would be 

y NP = |W F f (P nom ,K)W ; | p =J P (0) (19) 

In Ref. [16], the authors showed that the behavior and flatness of the effectiveness curve versus the estimation error is a 
good qualitative concept of robustness against that error; and if quantified, yields a good measure of robustness. The same 
reasoning, if applied to / P (<S), can lead to the following index J R as a suitable index of robustness: 

\ 1/2 
-^j-Et/p^WKO)] 2 (20) 

" d + j = o J 

where 5j= — l+2(j— l)/(n^— 1) for j=\ ...n s is a descritization of the uncertain parameter <S over its allowable range. 

Due to the conflicting behaviors of robustness and nominal performance, as has been shown in [16], a good strategy is 
to formulate a two-objective optimization problem. Based on the indices proposed in this paper, this would be 

|7np] 
mim } s.t. design constraints (21) 

x I Jr J 

where the solution for x (vector of design variables) will not be a single point x* only; but a set of Pareto optimal points. A 
point x* in the feasible design space S is called Pareto optimal if there is no other point x in the set S that dominates x* 
(improves at least one objective function without degrading another one) [16]. An extremely appealing characteristic of GA 
is the potential to make the population converge to the Pareto optimal set incorporating a great number of possible non- 
dominated compromises between the two objectives. This is a further motivation, in addition to the reasons discussed in 
Section 3, to adopt GA in the current investigation. Details on multi-objective GA concepts and implementation can be 
found in literature [40], as well as different examples of implementation, also including TMD system design [9,16,29]. 

Another important issue is the concept of worst-case performance [5,9] which motivates the formulation of the 
following index based onJp(<5) to characterize robust performance: 

7wcp= max J P (S)~ max J P (Sj) (22) 

6 allowable 1 <j < n$ 

The worst-case performance index in Eq. (22) guarantees avoidance of overestimation and thus achievement of high robust 
performance, if utilized as the single objective or as the preference criterion to select a solution from all those obtained 
through the two-objective optimization (21). As will be discussed in Section 5.1.1, the two mentioned ways of utilizing 
Jwcp are very similar in terms of the final result. Therefore, if robust design is the main purpose, single objective worst-case 
performance minimization is preferred over two-objective minimization (21) due to the much lower computational cost 
(single objective GA runs need much smaller population sizes and less number of generations). The other appealing 
characteristic of the worst-case performance, specifically in the case of p = oo, is that it can be generalized to the multiple- 
parameter case in a computationally efficient way. This is what we shall pursue in the next subsection. 

4.2. Multiple-parameter uncertainty 

Although the robustness concepts of the previous subsection can be directly generalized to the multiple-parameter 
uncertainty case, the exhaustive search over a fine descritized grid of the uncertain parameters hypercube is not a 
computationally acceptable method for more than two uncertain parameters. In this section, we shall focus on a 
computationally efficient way to assess the H x worst-case performance index for the case of generic multiple-parameter 
uncertainty. Here, the theoretical background is much more involved [14,30] and a more complicated uncertainty 
modeling is required. Similarly as we did in Section 2 to isolate the TMD system dynamics from the plant, now we need to 
pull out the uncertainties in the plant parameters and collect them into the structured perturbation block A within the LFT 
framework of Fig. 6. 

According to this LFT interconnection, the mapping from excitation v to output e (we assume that any input or output 
weighting is absorbed into P) is described as 

J A ev (s) = F U (F { (P,K),A) = ^ U (M,A) = M ev + M er \(l-M hr \)- x M hv (23) 
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Fig. 6. LFT framework with the TMD system as the controller and the structured uncertainty. 



Table 1 

Multiple-parameter uncertainty illustrative scenario. 



Uncertainty in the mass of all the floors 

Uncertainty in the stiffness of the first floor 
Uncertainty in the structural damping ratio 



m, — m,(l +Wj5j); i= 1 .2,. . .,N or in matrix form: 

M S =M S (I+W M A M ) 

where W M = diag(Wi,. . .,w N ), A M = diagfdi,. . .,S N ) 

/<, =k,(f+w N + 1 (5 N + ,) 

C = C(l+w N + 2 <5 N + 2) 



-1 <<5i<l for i=l,2. 



, JV+2; w,'s represent the uncertainty level; and the bar symbol represents the nominal value. 
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Fig. 7. Block diagram of the uncertain plant with uncertainties embedded in M s \ K s and £1. 

where a consistent partitioning of M=J r f (P,K) is assumed. JFJM.A) has a useful interpretation: M ev is the nominal mapping 
from v to e accomplished by the nominal structure-TMD system, and this nominal mapping is perturbed by A; M hr , M hv , and 
M e , reflect a priori knowledge on how the uncertainty affects the nominal map. Since the mapping from disturbance to target 
outputs includes all information required for characterizing the structure-TMD system performance, this LFT formulation 
seems quite promising. In particular, we are interested in calculating the following worst-case performance index: 



iv 



max 

allowable 



K(j(o)\\ 



max \\F u (M(jco),\)\\ 

A allowable" 



(24) 



At this point, two questions arise: 

1 ) How to cast a given problem into this LFT framework? 

2) How to calculate y W cp, >-. after having derived M and the structure of A for our system? 



The answer to the first question is a procedure referred to as pulling out the S's that we shall illustrate through a generic 
example, while the second is related to mathematical theorems and numerical algorithms in the field of robust control that 
we will briefly discuss and cite at the end of this subsection. 

Let us consider the N-story shear building with an n-MTMD on the top floor subjected to only base acceleration of 
unknown frequency content, i.e. W,= [l 0^ xN ] T , with only the inter-story drifts as target output quantities, see Eq. (15). 
Now assume the quite generic example of multiple parametric uncertainty described in Table 1 . 

It is easier to illustrate the procedure via block-diagram graphical representation. So consider the following 
rearrangement of Eq. (4) with a little change in notation, and its block-diagram representation in Fig. 7: 

ii = -aui-M,r 1 K s (/?cu+u)+M,r 1 G f F c -rgt' 

e = W u u, y = G a ii + v (25) 

where F C =F MTMD , v=a g , (i=2j{cox+co 2 ), and a=flaj- l m 2 . Note that G a and G F have been absorbed from Eq. (8) for the MTMD 
into the plant, so that F c and y are scalars and the MTMD block to be used as the controller K is SISO (G and G F should be 
omitted from matrices B K and C K of the MTMD block in Eq. (8)). 
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Fig. 8. Pulling out <5's from Mf 1 , K s , and fl blocks. 



Now consider the blocks that embed the uncertainties (highlighted in Fig. 7). It is an easy task to check that these blocks 
can be expanded as shown in Fig. 8 to separate the nominal values and the uncertain perturbations. For example, Fig. 8b is 
the graphical equivalent of the following relation: 



k,(l+Wjv + ic5 N+ i)+k2 
-k 2 





-k 2 
k 2 + k 3 







l<N 



= (I+W K2 <5 N+1 W K] ) 



k t +k 2 

-k 2 





-k 2 

k 2 +k 3 







(26) 



where W K1 = [1 1 ... 1] and W K2 = [w N+1 ... 0] T . Also note that in Fig. 8c, W~w N+2 \ N . 

By combining Figs. 7 and 8, and graphically rearranging the signals and blocks, the following block diagram is obtained. 

As can be seen, this block-diagram matches the LFT framework of Fig. 6. Accordingly, the structure of the uncertainty 
perturbation is A=blkdiag(diag((5 1 , ..., (5 N+1 ),(5 N+2 I N ). Also, the model of the uncertainty pulled out plant with state vector 
x=[u T ,ii T ] T , three groups of inputs r = [r M ,rJ,rT] T , v, and F c , and the three groups of outputs h=[hM, hj<, hJ] T , e, and y is as follows: 



(27) 



It is easy to verify that all the matrices without a subscript r or h are defined based on the nominal values as if there was no 
uncertainty, while those with these subscripts can be derived as follows: 
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Now we shall address the second question posed at the beginning of this subsection. Let us take a closer look at Eq. (24) 
for J WC p >00 and do some manipulations: 



Jwcp,oo = a max ll^uCMO'w^A)!^: 

A allowable 

= max max cf[jF u (M0'w),A)] : 

co A allow 



: max maxcr^utMO'co)^)] 

A allow co 



max 

1 < i < N u 



max a[F u (M(ja>i),A)] 

A allow 



(28) 



where a denotes the greatest singular value. First the maximizations have been interchanged, then the frequency range is 
discretized to N m grids. At each frequency cu,-, M(joj t ) is a constant complex matrix. This implies that we can focus on the constant 
matrix problem (the maximization inside the parentheses in the last expression), and re-solve at each frequency. In the robust 
control literature, there are investigations zooming on approximate solutions for the constant matrix problem. A well-known 
algorithm [30] uses the divide and conquer strategy based on upper and lower bounds for this problem. The lower bound is 
computed through an iterative procedure that mimics Hamiltonian methods for state-space norm calculations in its substeps. The 
upper bound calculation, on the other hand, is based on a LMI formulation. To obtain Jwcp.oo. such a procedure should be repeated 
over the frequency grid. However, since the peak-over-frequency is of interest, large frequency ranges can be quickly eliminated in 
order to reduce the computation time. Ref. [41 ] discusses relevant software implementation. 

The formulation developed in this section unifies the robust design of passive TMD systems for structures with multiple 
uncertainties with the rigorous SSV (fj.) framework in robust control (see Ref. [42] for a relevant example of active 
structural control). Besides real parametric uncertainty, this framework can also handle other types of uncertainties such 
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Fig. 10. FRF of the uncontrolled building from ground ace. to inter-story drifts of the first (a) and top floors (b). 

as the neglected dynamics introducing complex uncertainty blocks into A. Of course, this type of uncertainty is more 
relevant to active control where high-frequency neglected dynamics can have detrimental effects on system characteristics 
such as stability and noise immunity. The advantage of this unification is the possibility to adopt previously developed 
methods and inspire from ideas in the robust control literature for design of passive vibration absorbers as has been 
demonstrated in this section (Fig. 9). 



5. Implementation examples 



This section demonstrates how easily the method can be implemented on various problems with different assumptions, 
and illustrates the post-processing step of the proposed framework (Fig. 1). 

The structural model considered for all examples of this section is a 7-storied shear building with the following nominal 
structural properties: k=[2.5,2.5,2.5,2 2 2 2] x 10 7 (N/m), m=5.5 x 10 4 x [1,1, ...,1] (kg) and f s =0.01. The building is 
excited only by ground acceleration, and its frequency response (in terms of the inter-story drifts of the first and top floors) 
is depicted in Fig. 10. The structural mode peaks and frequencies have been also clearly illustrated. 

In the GA optimization, the non-dimensionalized parameters of the TMD units have been used as design variables. 
Assuming equal distribution of the total TMD system mass, the design variables vector x will comprise the frequency ratios 
Pi and damping ratios f j of the TMD units which can be conveniently searched within reasonably predictable bounds. 
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The parameters of the TMD units in the formulation of Section 2.2 can then be computed as 
m u = fr x m Slt or, kr.i = m T 0i(a S ] ) 2 , c u = 2£j^m Tii k TJ , fi t = fi tot /n 



(29) 



where n is the total number of TMD units; /i t ot is the ratio of the total mass of the TMD system to the total mass of the 
building m Sito t; and ctVi is the first natural frequency of the building. 

5.1. Single-parameter uncertainty examples 

Several design examples will be considered in this subsection, each having a different combination of TMD system 
specifications and minimization objectives. Of course, all the examples assume the single-parameter uncertainty (17) with 
W/ ; =0.05. All these cases (designs I, II, ..., V) are listed and summarized in Table 3 along with their respective assumptions. 
Each design example requires a two-objective GA run as described by Eq. (22). For instance, design I that is the typical 
central example (with which all other designs are compared) has been obtained by a GA run described in Table 2. 

The minimization objectives J NP and ] R are defined based ony P ((5) defined in Eq. (18) applying Eqs. (19) and (20), 
respectively. According to the designer's purpose, different types of norm, and input/output weights can be assumed for 
Jp(<S). In this study 3 different definitions for J P (<5) are considered as listed in Table 3. 

5.1.1. Different TMD system specifications 

To study the influence of the TMD system specifications on its effectiveness-robustness characteristics, in this 
subsection, an IMTMD on top floor (Design I), a UMTMD on top floor (Design II), and the "STMD on every floor" 
configuration (Design III) are designed as described in Table 3. The resulting Pareto fronts in the criterion space are 
compared in Fig. 1 1 . 

The Pareto front of an STMD on top floor has been also included in Fig. 1 1 for further comparison. As can be seen, the 
UMTMD system (with 3 design variables) dominates the STMD system (with 2 design variables), mostly in terms of 
nominal performance. Similarly, the IMTMD (with 14 design variables) dominates the UMTMD; this advantage is most 
significant in the high robustness (low J R ) region. In other words, irregular design enables us to design a highly robust 
MTMD system with considerably lower sacrifice of nominal performance. Finally, Design III dominates all previous designs 
by exploiting the "STMD on every floor" configuration that physically allows the installation of an overall heavier TMD 



Table 2 

Specifications of GA run for design I. 



Algorithm type 



Two branch tournament 



Minimization objectives 
Population size 
Pareto fraction 
Number of generations 
Cross-over 

Type 

Prob. 
Mutation probability 



7np andjK defined based on ]p\(5) 
500 

0.4 
400 

Scattered 

0.8 

0.05 



Table 3 

Specifications of design examples of Section 5.1. 



TMD system specifications 



Minimization objectives J m , and J R defined based on 



Irregular MTMD, n = 7, /( tot = 0.03 x= /3, ■•• ft fi 

Uniform MTMD, n = 7, ,Utot=0.03 x= \Afi ft jl ' 

STMD on every floor, /i tot =0.06x= [ft ■■■ ft Cl 

Irregular MTMD, n=7, /( tot = 0.03 X= [ft ... ft fi 

Irregular MTMD, n = 7,/( tot = 0.03x= [ft ■■■ ft fi 



Jn(S) = 


|W ,^ f (P,> 


.10 
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N!< i 


Jn(6) 








/pi(5) 








1P2(S) = 


|Woi^(P, 


.10 


Oivxi 


; P 3(<5) = 


[WoiJvfPj 


,10 


~w eq2 

Oivxi 



5.1.1/2 

5.1.1 
5.1.1 
5.1.2 

5.1.2 



a UMTMD— A0: total range of TMDs frequency ratio, O : offset of the central TMD frequency ratio /?, = 1 +/!„ -A/3/2, ft = /?,_, +A0/n(i = 2. . .7), 

C, = C(i=i...7). 

b Wot as defined in Eq. (15). 

c W eq , =W e[ ,|((a g = 3rad/s, C s = 0.3] (Eq. (13)). 

d W eq2 = lV eq |{a) g = 12rad/s, C g = 0.6) (Eq. (13)). 
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Fig. 11. Visualization of design sets in criterion space based onjp^d). 




Fig. 12. Behavior of Jpi(<5) against S for three design solutions. 



system on the building (a numerical value of ^ tot =0.06 compared to /i tot =0.03 is used in these examples). It is important to 
point out the significant effect of /i tot : if the "STMD on every floor" system were designed with /i t ot=0.03, it would be 
dominated even by the STMD on top floor. 

Now, let us analyze the trade-off between nominal performance and robustness by selecting three design points U, I 2 , 
and / 3 on the Pareto curve / as shown in Fig. 1 1 . These design points correspond to specific IMTMD system designs: /i is the 
most effective; J 3 is highly robust; and, I 2 achieves the best worst-case performance, i.e. the least value ofy WCP in Eq. (22) on 
the Pareto curve I. The/p(<5)-curves corresponding to these design points are illustrated in Fig. 12. 

As can be seen from Fig. 12, design point I lt being the most effective design under nominal conditions, is the most 
sensitive to uncertainty. On the other hand, the curve corresponding to design point J 3 is the flattest (the most robust in 
terms of J R ); but, it is less effective than J 2 in the entire uncertainty range of interest. In this sense, the design point I 2 is 
preferred to I 3 for achieving robust performance in the presence of uncertainty. Accordingly, minimizing the worst-case 
performance Jwcp is an appropriate preference criterion for choosing a robust solution from the entire set ofJ NP -J R optimal 
set. This would yield a solution close to the solution obtained by direct single objective minimization of/ WCP . Indeed, the 
interesting advantage of two-objective optimization is that it results in an entire set of Pareto optimal solutions that 
incorporate a broad spectrum of design options. Of course, the designer should decide whether or not obtaining such a 
spectrum of designs justifies the augmented computational cost of two-objective optimization. 

To further illustrate the characteristics of the two preferred TMD systems designs for achieving nominal (7i) and robust 
(J 2 ) performance, the uncertain frequency responses (from ground acceleration to the inter-story drifts of the first and top 
floors) have been depicted in Fig. 13. This figure shows that the robust design manifests lower sensitivity to the variation of 
the uncertain parameter (of course, the difference is not dramatic because the level of uncertainty is low). In order to find 
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Fig. 13. Uncertain FRFs (from ground ace. to the inter-story drifts of the first and top floors) for two designs It and I 2 (thick black curve: nominal response 
and red curves: uncertain samples). 



out the amount of vibration suppression accomplished by TMD systems 1^ 
Fig. 10 that shows the uncontrolled building FRF. 



and I 2 , it is useful to compare this figure to 



5.3.2. Different frequency contents of excitation 

A priori knowledge of the exact disturbance signal (or signals) is not available for design of the TMD system. Indeed, 
mitigation of the unknown disturbance is the primary task expected from the TMD system as a vibration controller. 
However, some information on the frequency content of the potential excitations might be available a priori; and we 
should try to take most possible advantage out of this information. For seismic vibration control, the ground excitation can 
be modeled by a stationary stochastic process with a spectral content of the form proposed by Kanai [43] and Tajimi [44], 
as previously presented in Eq. (13). In this subsection, the advantage of introducing such information into the design 
process will be investigated. Two spectra denoted by VV eq i (co g =3 rad/s, ( g =03) and W eq2 (to g =12 rad/s, £ g =0.6) in Table 3 
(with ground frequencies near the first and second structural frequencies) are considered in this study, as shown in Fig. 14. 
These spectra [33] match the frequency contents of Uemachi (a simulated ground motion using fault rupture model) and 
El-Centro (the well-known seismic event recorded during the Imperial Valley, California earthquake of May 18, 1940); and, 
are utilized in the definition of J P2 (S) (Designs IV) and J P3 (S) (Design V), respectively. 

The Design sets IV and V in/ NP -/ R criterion spaces based on7p2(<5) (with W eq i as the input weight) and J P3 (<S) (with W eq2 as 
the input weight) have been respectively shown in Fig. 15a and b, along with Design set I for comparison. At the first glance, 
these visualizations might seem ambiguous and require clarification. For each design set, two criterion spaces can be 
considered: one utilized for optimization, and the other used for visualization. Consider Design set IV, for instance, where J NP 
andy R based onJ P2 (S) have been utilized as objectives for optimization. In Fig. 15a, this design set has been visualized in the 
same criterion space as has been optimized; hence, it constitutes a Pareto front of non-dominated solutions in this space. In 
Fig. 15b, on the other hand, Design IV has been visualized in a different criterion space (the criterion space based on_/ P3 (<5)); 
consequently, it does not constitute a non-dominated Pareto front in this space. The opposite holds for Design set V; while, 
design I, in both subfigures, is visualized in criterion spaces different from the optimization space. Obviously, a post- 
processing (re-computation of objective indices) is required for visualization of a design curve in a criterion space different 
from its original optimization space. This post-processed visualization enables us to exactly assess the improvement of the 
effectiveness-robustness characteristics of the TMD system due to appropriate consideration of excitation spectral content 
and vice versa, i.e. the degradation effect due to not considering or inappropriately accounting for the frequency content of 
the disturbance. 
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Fig. 14. Kanai-Tajimi filters used in Designs IV and V. 
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Fig. 15. Visualization of design sets in criterion space based on (a)7p 2 (^) an d {b)J P3 (S). 



Compared to Design IV, Design V is closer to Design I in both criterion spaces (Fig. 15 a and b); in other words, influence 
of VV eql on the design is more significant than W eq2 as its shape shows more deviation from the constant (frequency- 
independent) filter W eq =\ used in Design I. 

Now, for final assessment and validation, let us analyze the efficacy of the designs by time-response simulations. 
The most effective points IV] and V] (from Designs IV and V) are selected and compared to the most effective design 
point I] (from the previous subsection) in terms of the first floor inter-story drift response under nominal conditions. 
For performance evaluation of these designs against seismic excitations with frequency content similar to the 
Uemachi scenario, response to stationary stochastic process has been simulated. Such a random signal is generated by 
filtering a white noise through W eq i and then scaling to achieve a peak ground acceleration (PGA) of 0.5 g. Fig. 16 
shows a sample of such a simulation. Final performance evaluation can be carried out by comparison of the peak and rms 
reduction of random vibrations that each design accomplishes. Table 4 presents the final results averaged over 1000 
simulations: 

As expected according to the previous observations from Fig. 1 5, the results of Table 4, by random vibration analysis in 
time domain, show that consideration of consistent frequency content (Design IV) improves the result compared to when 
no frequency information is taken into account (Design I) or inconsistent information is used (Design V); where in the 
latter case, the result is even worsened. This statement proves again true (by interchanging IV<->V) for the case of W eq2 
frequency content. This time validation is carried out via deterministic simulation against the El-Centro historical 
earthquake record, in Fig. 17. Although the improvements are moderate, the designs cannot be considered truly optimal if 
these improvements are not applied. 
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Fig. 16. A sample time-response to random noise colored by Uemachi filter, (a) Design point [, — peak reduction: 32.3%; rms reduction: 53.7%. (b) Design 
point IV, — peak reduction: 37.5%; rms reduction: 54.6%. (c) Design point V, — peak reduction: 30.4%; rms reduction: 53%. 

Table 4 

Random vibration analysis results (Uemachi). 



Reduction of: 


Design I 


iTO 


Design IV, (%) 


Design V, (%) 


Response peak 
Response rms 


35.7 
40.4 




37.1 
42.2 


34.6 
39.3 



(a) 

0.04 



-0.04 







Unctrl -I 


1 111 la 1 




-4\ 


\\^fif0myifi\ 


Vy^yA^^WV^ 



„(b) 

E 



Q 



0.04 



-0.04 

(C) 
0.04 



-0.04 




10 



15 20 25 

time [sec] 



30 



Unctrl 

IV, 






Unctrl , 

v, 




-^}f^ 


A»/ l \/V wv ^/\ 



35 



40 



Fig. 17. Time-response to El-centro earthquake record, (a) Design point U — peak reduction: 29.9%; rms reduction: 39%. (b) Design point IVi — peak 
reduction: 29.1%; rms reduction: 36.8%. (c) Design point Vt — peak reduction: 31.4%; rms reduction: 39.7%. 
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5.2. Multiple-parameter uncertainty example 

Consider the 7-storied shear building subject to the multiple-parameter uncertainty described in Table 1 with 
Wi=w 2 = ■■■ =w 7 =0.1, w g =0.2, and w 9 =0.75. Let us set up a single objective GA for designing an irregular 7-MTMD 
(14 design variables), once with/ NPo= as the single objective (nominal design), and once with_/ WCP ,. (robust design). All the 
assumptions and formulations follow the example of Section 4.2. The resulting MTMD designs (the damping and frequency 
ratios of the 7 TMD units) have been compared in Fig. 18. Due to the irregularity of the design and several number of 
influencing factors, specific results should not be generalized. However, there is a reasonable trend that can be generalized: 
more robust MTMD designs cover a broader frequency range and tend to have more damping. 

The FRF of the building equipped with each individual MTMD design has been shown in Fig. 19. The FRF of the robust design 
manifests lower sensitivity to uncertainties, especially near the second mode of vibration which is the critical frequency. 

Table 5 lists the numerical value of the nominal and worst-case performance indices for the two designs. As can be seen, 
the nominal design can suffer up to a 200% increase of gain (worst-case gain compared to the nominal gain) due to 
uncertainties; while, the gain increase for the robust design is at most 46%. Interestingly, for the building controlled by the 
robust MTMD, the critical frequency associated with the worst-case gain is near the second mode of vibration, while 
the nominal peak gain occurs at the first mode. For the building with the nominally optimal MTMD both peaks occur near 
the second mode. 

The worst-case combination of parameters happens for <5i=<5 2 =<55=(S7=l and <5 3 =<5 4 =(55=<5 8 =<5 9 = — 1. It corresponds to 
the least allowable stiffness of the first story and structural damping, and a distribution of the mass of the floors that favors 
the second vibrational mode-shape (least allowable mass for the 3 floors in the middle, and the maximum allowable for the 
others). An interesting fact about this worst-case scenario is that it has happened at a corner of the uncertain parameters 
cube. Based on experience with several numerical simulations, the authors suspect that this is true for any problem of the 
generic class considered in this paper. This can avoid the cost and complexity involved with computing the worst-case gain 
bounds, and would lead to the exact value. The worst-case combination of parameters that the existing algorithms, e.g. the 
function wcgain in MATLAB Robust Control Toolbox [45], report is just the parameter combination attained at the end of 
the iterative search procedure for computing the lower bound on the worst-case gain. For the problem at hand, this will 
not always yield a corner; but, numerical experiments on this problem have shown that searching only the corners always 
improves (or at least, maintains) the result obtained from the lower bound procedure. However, no relevant theoretical 
grounding has been found in the literature. Hence, if theoretical guarantee is required, the divide and conquer strategy 
based on upper and lower bounds must be carried out. 
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Table 5 

Comparison of the nominal and robust designs in terms of performance indices. 
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Fig. 19. Uncertain FRFs (from ground ace. to inter-story drifts of the first and top floors) for nominal and robust designs (thick black curve: nominal 
response and red curves: uncertain samples). 

6. Conclusions 

Toward establishing the proposed generalized framework for robust optimal design of TMD systems, the following 
primary contributions and results have been presented: 



1 ) The LFT formulation: A rigorous LFT formulation has been introduced that isolates the TMD system as a feedback controller 
and the uncertainties as a structured block perturbing the nominal dynamics. It has been shown that this formulation can 
accommodate a general structure with any TMD system configuration and multiple-parameter uncertainty. 

2) Generalized performance indices for optimization: The concept of weighted MIMO norms from external disturbance 
inputs to control-target quantities has been used as the most general measure of performance. In this way, knowledge 
about the frequency content of the excitations can be incorporated and the H 2 or H&, optimization of structural 
response in terms of different quantities (inter-story drifts or floor accelerations) can be formulated in a unified manner. 
A nominal performance index, a robustness index and a worst-case performance index have been formulated 
accordingly. The latter index in conjunction with the structured uncertainty modeling strategy has been shown to be 
efficiently computable in the presence of multiple-parameter uncertainty. 

3) Exploitation of evolutionary heuristics: Advantages of exploiting evolutionary heuristics were demonstrated through 
implementation of GA optimization. 

4) Assumptions on the TMD system: Comparison of IMTMD and UMTMD solution sets, designed to be installed on the top 
floor of a multi-story building, have demonstrated domination (in terms of effectiveness-robustness Pareto optimality) 
of irregular design, quite significant in the high robustness region. It has been also shown that by spatial distribution of 
the TMD units using a simple configuration of "STMD on every floor", a totally heavier and hence dominant TMD system 
can be exploited. 

5) Frequency content of excitation: Consideration of frequency content information, through application of input weighting 
filters, has been shown to improve effectiveness-robustness characteristics of the TMD systems against excitations of 
consistent spectral content. Accordingly, a seismological study on the site plan of implementation and an estimation of 
the wind loadings spectral content can be beneficial to the TMD system design. 

6) The active control perspective toward passive optimization: Application of ideas in active control theory for passive design of 
TMD systems has proved particularly beneficial: it has enabled the development of a flexible and rigorous formulation 
and utilization of worst-case performance assessment algorithms previously developed in robust control literature (far 
more superior than exhaustive search and existing Monte Carlo methods available for passive optimization). In a more 
general sense, the present paper can be considered a first step toward unification of passive and active TMD systems 
design in view of robust control. Future research in this context can involve simultaneous synthesis of the passive 
mechanical parameters and the control law for active TMD systems, or utilization of fixed-structure active control 
synthesis schemes for design of passive TMD systems (e.g. instead of the GA optimization scheme). 
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Appendix. Formal definitions and computation of the weighted MIMO norms 

Consider the general input-output weighted LFT configuration of Fig. 5 formulated as the cascade interconnection of 
the input weight, the non-weighted LFT and the output weight: 
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according to the state-space realizations of the input and output weighting filters; the state-space matrices for the (non- 
weighted) LFT configuration described by Eq. (3-II); and, the resulting augmented matrices for the weighted LFT 
configuration, described as 
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This augmented system T ev (s), a stable structural plant controlled by a passive mechanical controller and weighted 
through proper stable filters, is a proper and real rational stable transfer matrix (s TIH^); hence, its H,, norm exists and is 
defined as [14]: 
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Re(s) > 60eR 

with a denoting the greatest singular value. Furthermore, T ev (s) is also strictly proper (e IZH2X provided that 

D fl = D D CL D ; = 
holds, allowing the following definition [14] for its (finite) H 2 norm: 
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Condition (A-3) is quite easily satisfied: considering the definition for D C l in (3-II) and substituting D zd from (5) and 
D K =0 for both TMD configurations (see Eqs. (8) and (9)) leads to 
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where we have used a block partitioning of numeric matrices D and D, consistent with the partitioning of the transfer 
matrices VV and W/ presented in Eqs. (14) and (11), respectively. Eq. (A-5) shows that (A-3) holds iff Do.accMs^D/ext is zero. 
This means that if the input-output path from external force disturbance to acceleration target output is weighted through 
an input or output strictly proper filter, (A-3) will be satisfied. Under this condition, the 2-norm is easily computed based 
on the state-space matrices after solving a Lyapunov equation [14] for which efficient algorithms exist (e.g. see Ref. [46]). 
The computation of the H x norm (not requiring Eq. (A-3) to be satisfied), on the other hand, is more involved and needs 
numerical search. Direct application of the formal definition equation (A-2) is not efficient; however, there are other well- 
established numerical algorithms, e.g. Ref. [47] presents an efficient algorithm (with guaranteed accuracy) based on the 
relation between the singular values of the transfer function matrix and the eigenvalues of a related Hamiltonian matrix 
constructed from the state-space matrices of the system. 
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